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ABSTRACT 

Observations of quasi-periodic oscillations (QPOs) in X-ray binaries hold a key to understanding 
many aspects of these enigmatic systems. Complex appearance of the Fourier phase lags related 
to QPOs is one of the most puzzling observational effects in accreting black holes. In this Letter 
we show that QPO properties, including phase lags, can be explained in a framework of a simple 
scenario, where the oscillating media provides a feedback on the emerging spectrum. We demonstrate 
that the QPO waveform is presented by the product of a perturbation and a time delayed response 
factors, where the response is energy dependent. The essential property of this effect is its non-linear 
and multiplicative nature. Our multiplicative reverberation model successfully describes the QPO 
components in energy dependent power spectra as well as the appearance of the phase lags between 
signal in different energy bands. We apply our model to QPOs observed by Rossi X-ray Timing 
Explorer in BH candidate XTE J1550-564. We briefly discuss the implications of the observed energy 
dependence of the QPO reverberation times and amplitudes to the nature of the power law spectral 
component and its variability. 

Subject headings: accretion, accretion disks — black hole physics — X-rays: binaries — stars: individual 
(XTE J1550-564) 



1. INTRODUCTION 

Galactic X-ray Black Hole (BH) candidates exhibit 
spectacular evolution of timing and spectral behavior 
(IRemillard fc McClintockl I2006T IBellonil 120051) . Strong 
variability is usually observed in states, dominated by 
non-thermal emission (i.e. hard and intermediate states). 
The variability properties st rongly correlate with over- 
all sp ectral state evolution (jKlein-Wolt fe van der Kiisl 
2008). A common property of the variability is nearly 
periodic flux modulations with frequency between 0.01 
and 20 Hz generally referred to as low frequency quasi- 
periodic oscillations (QPOs). The exact origin of QPOs 
in BH systems is still under debate. In this Paper we 
present a new important evidence which allows us to con- 
strain the QPO models and to advance in understanding 
of this phenomenon and as well as the overall BH phe- 
nomelogy. 

QPOs are observed as one or more relatively narrow 
peaks in Fourier power density spectra (PDS). QPO sig- 
nal is shown to be strongly c oupled with the power law 
part of the energy spectrum (jShaposhnikov et all 120111 : 
ISobolewska fc Zvcki 1 120061 ). Probably the most intrigu- 
ing aspect of the QPO phenomenology is the behavior 
of Fourier phase lags, i.e. differences in phases calcu- 
lated in two different energy ranges. High time resolu- 
tion data on brightest BH transients, provided by Rossi 
X-ray Timing Explorer {RXTE) allowed to study this 
effect in detail. Phase lags observed during the ini- 
tial intermediate state of 1998 event o bserved in XTE 
J1550-564 have shown complex behavior (|Cui et al.ll2000l : 
IRemillard et al.l I2002D . Particularly remarkable is the 
fact that the fundamental QPO and its first harmonic 
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exhibit lags of opposite signs. While the fundamental has 
shown the negative (or soft) lags, the lag in the first har- 
monics was hard for most observations (conventionally 
we define positive or hard lags when signal is delayed at 
harder energies and vise versa for negative or soft lags). 
Similar behavior was observed during 1999 outburst from 
XTE J1859+226 (jCasella et al.ll2004D a nd in persistent 
BH candidate GRS 1915+105 (jLin_et^aI|[2003). The con- 
clusion of these studies was that either the lags behavior 
is due to different physical origin of QPO harmonics or 
some unknown physical mechanism lead to a particu- 
lar QPO waveform dependence on energy, which is well 
reproduced in many observations. However, no QPO 
model have met the challenge to account the observed 
phase lag behavior. 

In this Letter we present new and extremely impor- 
tant paradigm which allows much better understand- 
ing of QPO behavior and put a stronger constrains on 
its physical origin. First, we present a QPO wave- 
form parametrization which provide proper description 
of QPO properties both in term of Fourier amplitudes 
and phases. In fact, we develop the first analytical model 
for QPO phase lags which allows direct modeling of the 
phase lags seen in observations. Secondly, based on the 
results of our model, we draw conclusions on the physi- 
cal nature of the observed oscillations and the oscillating 
media. We tentatively call our model the "QPO rever- 
beration model" . We do not introduce any new physical 
paradigm. Instead, acting on the premises of the stan- 
dard truncated disk/hot inner corona model, we identify 
a specific way, in which all major observed QPO prop- 
erties are naturally produced. The main idea behind 
the "QPO reverberation" model is that the inner part 
of the accretion flow acts as a non-linear forced oscilla- 
tor where both perturbation and response signals enter 
the observed signal in multiplicative manner. While the 
perturbation is hydrodynamical in nature, the response 
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factor reflects the spectral shape reaction to the changes 
in oscillating media properties introduced by the initial 
perturbation. 

The QPO reverberation model matches the data very 
well, allowing to convert a complicated form of the phase 
lags into well understood physical delay time. The re- 
sponse signal is always delayed with respect to the per- 
turbation, as expected in case of forced oscillations. The 
main conclusion which we draw based on our analysis 
is that multiplicative reverberation is a primary mech- 
anism leading to the non-sinusoidal nature of the QPO 
waveform and defining the behavior of QPO phase lags. 

In the next section we describe the model and our 
RXTE data analysis. In §3 we briefly discuss the impli- 
cations of our results to different QPO scenarios. Con- 
clusions follow in the §4. 

2. REVERBERATION EFFECT IN QPO 
2.1. Model 

We parametrize the QPO signal by the following wave- 
form: 

s(t) cx (l + ae~ A * cosw i)(l + &e _At cos(w (t-t d ))). (1) 

We assume that an oscillation process in the system 
signified by QPO leads to a time delayed response of the 
spectral form of the outgoing radiation. The motiva- 
tion for this functional form of the signal is most easily 
understood In terms of the power law shape of the av- 
erage energy spectrum NE~ a . If we assume that the 
normalization N and the index a experience exponen- 
tially decaying oscillations, where oscillations in index 
are delayed with respect to the oscillations in the nor- 
malization then, after expanding the power law in Tay- 
lor series keeping only the first term, we arrive to the 
expression (fT]). The Fourier transform of the waveform 
and its products (i.e. power spectrum, phase lags) are 
presented in Appendix. 

2.2. Application to data 

We apply our model to the data collected with RXTE 
during the brighte st BH transient event observed from 
XTE J1550-564 dHoman et alJ l205lt ICui et al.l |2(M 
Sha poshnikov fc Titarchukl 12009ft . We use a combina- 
tion of Binned and Event data modes to extract 7.8125 
millisecond time resolution lightcurves in channels cov- 
ering 1.5-3.5, 3.5-5.0,5.0-6.0, 6.0-7.5, 7.5-12.0, 12.0-14.6, 
14.6-19.3 and 19.3 - 50.0 keV energy ranges. We com- 
pute Fourier transforms for consecutive intervals of 64.0 
sec for each lightcurve. We then compute correspond- 
ing power spectra and phase lags as a complex argument 
of the cross spectrum using the first channel range as a 
reference. We rebin these data products using 1.02 log- 
arithmic rebinning factor. In Figures Q] and [2] we show 
representative power spectra and phase lags for two ob- 
servations, 30188-06-05-00 (Observation 1) and 30191- 
01-30-00 (Observation 2). The fundamental QPO during 
these observations appeared at 2.35 Hz and at 6.5 Hz re- 
spectively. The first observation is taken during a harder 
spectral state, i.e. earlier in the outburst. It showed pos- 
itive lags for both QPOs, while the second observation 
shows negative lag at the fundamental. For each ob- 
servation we combine eight PDSs and seven correspond- 
ing phase lags spectra in a joint model fit. Using the 
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Fig. 1 .— RXTE data from observation of XTE J1550-564 (Ob- 
sID 30188-06-05-00). The bottom panel shows phase lags between 
signals in soft and hard energy ranges. The middle panels show 
PDSs for the two energy ranges. The data are fitted with the re- 
verberation model, shown as solid color lines (red for the soft range 
and blue for the hard). For the PDS fits the dashed lines show con- 
tribution of Poisson and broad-band noise components. The top 
panel shows the QPO waveforms for signal for both energies the 
soft (green) and the hard energy ranges as given by the model fit 
parameters. 

QPO waveform described by Eq. [TJwe model QPOs in 
the power spectra and corresponding phase lags between 
signals. The parameters describing perturbation signal a 
and A are the same for all waveforms, while the parame- 
ters describing spectrum reverberation, i.e. its amplitude 
b and time delay td (or phase delay <p — 27rajn^d) are in- 
dividual for each energy range. The data and the model 
fits for two energy ranges are presented in Fig. [T] and [21 
Some systematic is seen in the phase lags fits in the 
vicinity of the first harmonic of the QPO in power spec- 
tra. It would be naive to expect that the parametriza- 
tion using the expression (fT]) will describe every aspect 
of the data. Rather than the exact QPO waveform, the 
presented model provides a basic framework for develop- 
ment of the future models which capture such effects as 
contribution of aperiodic PDS component in the phase 
lags, effects of statistical noise, etc. Taking this fact, 
the model fits the data extremely well capturing major 
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Fig. 2.— The same as Fig Q] for RXTE observation 30191-01- 
30-00 which is performed later during the outburst. Note that for 
this case, as opposed to the case presented in Fig[T]the phase lags 
at the fundamental and the first harmonic have different signs. 

behavior of the phase lags and the exact form of the 
fundamental QPO, which is energy dependent and assy- 
metrical. In Figures [3] and U] we show the dependence of 
the reverberation times and amplitudes b on energy. 
Below we discuss the implications of these results for the 
physical picture of QPO production in BH candidates. 

3. DISCUSSION 

The results of the presented model fits to the data show 
that the reverberation scenario properly describes the 
QPO properties. Apparent non-linear multiplicative na- 
ture of the effect and its energy dependence naturally ex- 
plains a number of observational effects including shape 
of the QPOs in PDS, energy dependence of the QPO 
powers, evolution and appearance of the phase lags. 

Probably even more important result of our analysis 
is the fact that negative Fourier time lags observed in 
many cases for the fundamental QPO are of a purely 
mathematical nature and should not be treated as physi- 
cal times. The reverberation model properly convolves a 
complicated behavior of Fourier phases into a one, well 
understood, physical time td which can then be consid- 
ered in terms of physical models. 
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Fig. 3. — Energy dependence of the reverberation time td for 
the Observation 1 (blue points, a = 1.30 ± 0.06, A = 0.56 ± 0.01) 
and Observation 2 (red points, a = 1.0 ± 0.1, A = 2.07 ± 0.03). 
The logarithm f(E) = A log(E) + B fit is shown by black solid line 
for the first observation and dashed line for the second one. The 
time delays for the first observations clearly follow the logarithm 
dependence while the data for the Observation 2 is consistent with 
the constant delays. 

Another important result is the nature of the ob- 
served reverberation and energy dependence of the re- 
verberation time and amplitude . It is clear that the 
perturbation actively changes the condition in the sys- 
tem. Geometrical QPO models which are based on a 
mechanical rotation or precession listeria fc VietrilHoM 
ISchnittman et al.1 120061 . see e.g.) lack this active ingre- 
dient. Reverberation effect strongly favors QPO mod- 
els which involve oscillations affecting the state of mat- 
ter which then is reflected in the form of the emerging 
spectrum. This grou p of models includes tr ansition layer 
(corona) oscillations (jTitarchuk et al.lll998f) and gravita- 
tiona lly trapped disk oscillation modes (Wagoner et al. 
|2001[ ). Logarithmic dependence of the time delays for the 
hard state (see Fig. [3] for Observation 1) strongly favors 
accretion induced oscillations of the Compton corona as 
model for QPOs. 

In the truncated disk scenario a thin accretion disk 
extends down to some inner radius where it superseded 
by a optically-t hin geometrically thi ck configuration or 
corona (see e.g. iTomsick et aTl 12009). In this scenario, 
the thin disk is responsible for the soft thermal part of 
the spectrum, while the corona produces a non-thermal 
power law part of the energy spectrum. Variation in the 
inner disk radius, caused by surges of accreting matter, 
introduce perturbation in the corona and modulate its 
physical parameters, i.e. its density, temperature, optical 
depth, etc. This leads to temporal variations of the spec- 
tral shape of the escaping radiation. The time delay be- 
tween the perturbation and the spectral response is set by 
sound speed in the corona. The logarithmic dependence 
of time delays with energy observed for Observation 1 
(hard state) indicates (see Figure [3]) , according to ex- 
pectations, that Comptonization play significant part in 
spectral formation during these states. Namely, multiple 
upscattering events experienced by photons cause wave- 
like structures to propagate from the low energies toward 
the higher part of the time dependent spectrum. On the 
other hand, during the soft intermediate state (Obser- 
vation 2), the reverberation time is flat with respect to 
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Fig. 4. — Energy dependence of the reverberation amplitude b for 
the Observation 1 (blue) and Observation 2 (red). The Observation 
2 data is fitted with f(E) = Alog(l/ E)+B which is shown by black 
solid line. 

energy, indicating that the spectrum reacting to the per- 
turbation at the same time at all energies. Moreover, the 
reverberation amplitude for this observation is well rep- 
resented by b oc log(l/.E), in agreement with variability 
being produced by a pure pivoting power law. This can 
not be facilitated through the thermal Comptonization 
and will require cither synchrotron or non-thermal (such 
as bulk motion, see e.g. Titarchuk&Shaposhnikov 2010) 
Comptonization mechanisms to be invoked. This subject 
is beyond this Letter and will be addressed in more detail 
i n upcoming public a tions. 

iMaccarone et al.l (|2011l ) studied QPOs observed in 
GRS 1915+105 by means of bispectrum and found evi- 
dence for strong non-linear coupling between QPO har- 
monics and the broadband noise and considered model 
involving non-linear oscillator. Our findings also strongly 



suggest that the inner accretion flow region acts as a non- 
linear forced oscillator, which, by means of combination 
of radiation mechanisms such as thermal/non-thermal 
Comptonization and synchrotron radiation, modulate 
the emerging spectrum of the system, leading to a non- 
linear nature of the QPO signal. The waveform presented 
by expression ([T]) predicts that on time scales shorter 
than the exponential decay time 1/A both the average 
flux and variability are proportional to e~ xt and, there- 
fore, are linearly related to each other. Rms-flux relation, 
well established in broad variability from BH sources, 
was recen tl y fou nd in QPOs from XTE J1550-564 by 
iHeil et al] ( 20111 ). A specific feature of the relation is 
that it sw i tches from positive to negative above 5 Hz. 
IHeil et alJ (|2011[ ) show that scaling of the QPO power to 
the power of underlying broad band noise corresponding 
to the QPO frequency reinstates the linear rms-flux re- 
lation, in excellent agreement with the expectation that 
the noise variability is associated with the driving force 
for non-linear oscillator associated with QPO signal. 

4. CONCLUSIONS 

We present a new observational effect in QPOs. We 
identify the time delayed reverberation in QPO signal as 
mechanism facilitating the phase lag behavior observed 
in bright galactic BH candidates. The QPO reverbera- 
tion model successfully describes the QPOs in energy de- 
pendent power spectra and the corresponding phase lags. 
The presented results strongly point to oscillation of the 
Comptonizing and, most probably magnetized, corona as 
the origin of the QPO signal from accreting BHs. The 
QPO reverberation parametrization presents, for the first 
time, a consistent description of the QPO signal both in 
Fourier amplitude and phase domains, providing an ex- 
citing avenue to finally resolving the QPO nature. 



APPENDIX 
QPO WAVEFORM 
Let us consider the waveform described by the expression: 

s(t,E) oc (1 + ae~ At cosw t)(l + b(E)e~ xt cos(cj < - <t>(E))), 

where b(E) and 4>(E) are functions of energy E. 
The Fourier transform of s(t) is given by 



S(u,E) = 



1 



a + be 1 



a + be 



abe z 



abe 



ab cos cf> 

2A + iui ' X + i(uj-uj ) ' X + i(uj + uj ) ' 2(2A + i(u> - 2cj )) ' 2(2A + i(u + 2cj )) 



(Al) 



(A2) 



The QPO power spectrum and phase lags between signals at energies E\ and E-i are calculated in straightforward 
manner, i.e.: 

|2 



P(E) = \S(E lt2 )\' 



and 

correspondingly. 



5P(E 1 ,E 2 )=Aig{S(u,E 1 )*S(u,E 2 )) 



(A3) 
(A4) 
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